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Abstract 

We calculate the decay rate of a Yukawa fermion in a thermal bath 
using finite temperature cutting rules and effective Green's functions ac- 
cording to the hard thermal loop resummation technique. We apply this 
result to the decay of a heavy Majorana neutrino in leptogenesis. Com- 
pared to the usual approach where thermal masses are inserted into the 
kinematics of final states, we find that deviations arise through two dif- 
ferent leptonic dispersion relations. The decay rate differs from the usual 
approach by more than one order of magnitude in the temperature range 
which is interesting for the weak washout regime. We discuss how to 
arrive at consistent finite temperature treatments of leptogenesis. 



1 Introduction 

Leptogenesis [112] is an extremely successful theory in explaining the baryon 
asymmetry of the universe by adding three heavy right-handed neutrinos TV; to 
the standard model, 

SC = iNid^Ni - X^iaNi^ia - ^MiNiN? + h.c. . (1) 

with masses Mj at the scale of grand unified theories (GUTs) and Yukawa 
couplings X Ui i a similar to the other fermions. This also solves the problem of 
the light neutrino masses via the see-saw mechanism without fine-tuning [3H5]. 

The heavy neutrinos decay into lepton and Higgs boson after inflation, the 
decay is out of equilibrium since there are no gauge couplings to the stan- 
dard model. If the CP asymmetry in the Yukawa couplings is large enough, 
a lepton asymmetry is created by the decays which is then partially converted 
into a baryon asymmetry by sphaleron processes. As temperatures are high, 
interaction rates and the CP asymmetry need to be calculated using thermal 
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field theory [BHH] rather than vacuum quantum field theory. However, in the 
conventional approach [7], thermal masses have been put in by hand without 
investigating the validity of this approach in detail. We have adressed this is- 
sue in [H| and found that corrections arise through the occurence of two lepton 
dispersion relations in the thermal bath. In this paper, we calculate the decay 
rate of the heavy neutrino in a consistent way (chapter [3]) which automatically 
includes the effect of leptonic quasiparticles, compare it to the conventional 
approach (chapter [5]) and give an outlook of what needs to be done to arrive 
at consistent descriptions of leptogenesis. Our calculation is general enough to 
be applied to all decays of a Yukawa fermion at finite temperature, which has 
interesting implications for other early universe dynamics (chapter [5]). 

2 Hard thermal loops and thermal masses 

If a particle reaction like scattering or decay takes place in the background of 
a heat bath, e.g. in the hot state of the early Universe, thermal field theory 
has to be employed to describe this process. There are two different approaches 
for considering finite temperatures within quantum field theory, the imaginary 
and real time formalism IIP] , both yielding the same results. In this work, we 
will use the imaginary time formalism. Going from zero to finite temperature, 
ensemble-weighted expectation values of operators have to be used rather than 
vacuum expectation values. For an operator A, this reads 

(A) p =tr(pA), (2) 

where p is the density operator describing the ensemble. In this way it can be 
shown that the propagator at finite temperature T is given by its usual vacuum 
expression where the zero component of the momentum is replaced by imaginary 
discrete Matsubara frequencies go = 2nmT in the case of bosons or (2n + l)i7rT 
in the case of fermions with integers n (see e.g. Perturbation theory at 

finite temperature then follows from using these propagators and summing over 
the Matsubara frequencies in loop diagrams. 

However, using these bare thermal propagators can lead to inconsistent re- 
sults, which are not complete to leading order, infrared divergent, and gauge 
dependent in the case of gauge theories. A famous example is the damping 
rate of a plasma wave in the quark gluon plasma, which is different in differ- 
ent gauges. In order to cure this, the hard thermal loop (HTL) rcsummation 
has been invented |12ll3j . For this purpose one has to distinguish between 
hard momenta of the order T or larger and soft momenta of the order gT or 
smaller, where g is the coupling constant, which is strictly possible only in the 
weak coupling limit j < 1. After all, the HTL improved perturbation theory 
has been successfully applied to thermal QCD for the description of the quark- 
gluon plasma (see e.g. [H]). The basic idea is that the bare propagators are 
replaced by resummed propagators, if the external momentum is soft Q < gT. 



2 



— ©— 



n W n 



Figure 1: Resummed propagator 



For a scalar field (Fig. [I]), this resummation follows from the Dyson- Schwinger 
equation as 

.A* = *A + iA(-m)iA + ■ ■ ■ = = Q2 _^_ n - (3) 

The thermal self-energy II of the scalar held then acts as a thermal mass m% h = H 
and gives a correction to the zero-temperature mass m% ot = tUq + m^ h - Since 
II is of the order ~ gT, the resummation will only affect the propagator when 
Q < II ~ gT, which is reflected in the prescription to resum only soft momenta. 
The resummed fermion propagator has a more complicated structure and will be 
explained in the next section. In general, the self-energy is momentum depen- 
dent, e.g. the photon self-energy in QED. In this case, the leading order gauge 
independent self-energy follows from integrating only over hard momenta in the 
loop diagram defining the self-energy. This HTL contribution in the resummed 
propagator leads to a correction of the order gT which cannot be neglected if 
the momentum of the propagator is soft. The poles of the HTL resummed prop- 
agators then describe the dispersion relations in the medium, e.g. plasma waves 
following from the resummed photon propagator. In addition to propagators 
also HTL effective vertices related to the propagators by Ward identities might 
have to be used. 



3 Decay and inverse decay rate 

In the neutrino decay we want to calculate, the Higgs boson and the lepton 
acquire thermal masses of the order ~ 0.2 — 0.4 T via their interactions 
with other standard model particles. In the regime where the temperature is of 
the order of the neutrino mass T ~ M, one of the momenta of the decay products 
can be soft and has to be resummed. In the regime where M < 0.2 — 0.4 T, 
both Higgs boson and lepton momentum will be soft and need to be resummed. 
We are interested in both regimes, therefore we will resum both Higgs boson 
and lepton propagator. The case of resumming only one propagator is included 
in this approach, since resumming a hard propagator gives only a negligible 
correction to the bare propagator. The HTL resummation has been invented 
for the weak coupling limit g <C 1. This limit does not apply in our case, our 
more phenomenological approach is rather motivated by the desire to capture 
effects beyond perturbation theory and justified a posteriori by the sizeable 
corrections it reveals, similar to the treatment of meson correlation functions 
in p3]. 
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Figure 2: N decay via the optical theorem with dressed propagators denoted 
by a blob 



Wc consider a leptogencsis-inspired model with a massive Majorana fcrmion 
JV coupling to a massless Dirac fcrmion t and a massless scalar <f>. The interaction 
and mass part of the Lagrangian then reads 

A„t,ma SS = gN4>i - ^MNN C + h.c. , (4) 

The HTL resummation technique has been considered in [16] for the case of a 
Dirac fermion with Yukawa coupling, from which the HTL resummed propa- 
gators for the Lagrangian in Eq. ^ follow directly. We like to calculate the 
interaction rate L of N <-> £<fr. 

We cut the N self energy and use the HTL resummation for the fermion and 
scalar propagators (Fig. [2]). 

According to finite-temperature cutting rules |17|18) . the interaction rate 
reads 

r(P) = -J-tr[(^ + M)Im£(P)]. (5) 
At finite temperature, the self-energy reads 

P L S*(K)P R D*(Q), (6) 

where Pl and Pr are the projection operators on left- and right-handed states, 
Q = P — K and we have summed over Majorana and Dirac spins. 
The HTL-resummed scalar propagator is 

where = g 2 T 2 /12 is the thermal mass of the scalar, created by the interaction 
with fcrmions. Due to the reduced Majorana degrees of freedom, differs from 
the Dirac-Dirac case by a factor 1/2 [To] . 

The effective fermion propagator in the helicity-eigenstate representation is 
given by [T9H2T] 

S*(K) = ±A + (K)( 7o - k • 7 ) + \&-{K){ lQ + k • 7 ), (8) 
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where 

A ± (K) = 

and 



±fc - k k + k 
- k0±k+ —{ ±l -^2k— lTl k^k 



(9) 



ml = ^T 2 . (10) 

This again differs from the Dirac case by a factor 1/2 [IB] . 
The trace can be evaluated as 

ti[(f> + M)P L S*(K)P R } = A + (p -pr)) + A_(po+pri), (11) 

where r\ = cos 6* is the angle between p and k. We evaluate the sum over Mat- 
subara frequencies by using the Saclay method [22] . For the scalar propagator, 
the Saclay representation reads 

f P j 

D*(Q) = - dre* T — {[1 + n B {uj q )]e-^ + n B {w q )e u " r }, (12) 

JO 

where f3 = 1/T, n B (^ q ) = Xjie^^ — 1) is the Bose-Einstein distribution and 
UJ q =c J 2 + m \- F° r the fermion propagator it is convenient to use the spectral 
representation [2"3"] 

A±(K) = - dr'e fcoT / duj p ± (u, k)[l - n F (u)]e- UT , (13) 

JO J-oo 

where np(u;) = l/(e LJ ' 9 + 1) is the Fermi-Dirac distribution and p± the spectral 
density [1^0] . 

The fermion propagator in Eq. ((8]) has two different poles for 1/A± = 0, 
which correspond to two leptonic quasiparticles with a positive (A+) or negative 
(A_) ratio of helicity over chirahty [24 27 . The spectral density p± has two 
contributions, one from the poles and one discontinuous part. Since the quasi- 
particles are our final states, we will set K such that l/A±(K) = 0. 

Thus we are only interested in the pole contribution 

P P ± lc (u, k) = _ W±) + S{UJ + W= ^' (14) 

where w± are the dispersion relations for the two quasiparticles, i.e. the solu- 
tions for fco such that 1/A ± (w±,k) = 0, shown in Fig. [3J There exists an 
analytical solution for u>± making use of the Lambert W function which has 
not yet been reported in the literature. The analytical solution is explained in 
detail in the appendix. One can assign a momentum-dependent thermal mass 
m±(k) 2 = uj±(k) 2 — k 2 to the two modes as shown in Fig. 2] and for very large 
momenta the heavy mode m+ approaches y/2me, while the light mode becomes 
massless. 

In order to execute the sum over Matsubara frequencies, we write fco = iu n 
with u! n = (2n + 1)ttT and remember that when evaluating frequency sums, 
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Figure 3: The two leptonic dispersion relations compared with the standard 
dispersion relation to 2 = k 2 + mj in blue are shown. 



also po — ito m = i(2m + 1)tvT can be written as a Matsubara frequency and 
later on be continued analytically to real values of po [10128129] , In particular 
-1. We can write 



e PO0 — gi^mP 



then 



n n / — — oo 

T ^ e (po-feo)r e fe r' =e P0T^ r /_ T ) j 



(15) 
(16) 



since — /3 < r' — r < j3. After evaluating the sum over fco and carrying out the 
integrations over r and r', we get 



/OO i 
dupiCw.fc) — 
2w 9 



1 + n B (uJ q ) - n F (io) 
pa - lo - Lo q 

| n B (u q ) + n F (u) 

PO - LO + LJ q 



(17) 
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Figure 4: The momentum-dependent quasiparticle masses m\ = w± — k 2 are 
shown. 



Integrating w over the pole part of p± in Eq. (fT4| , we get 



feo 



2w„ 



-k 2 



2m 2 
2 k 2 



1 + tib — n F 



n>B + n F 



2mi 



p -UJ±- LO q po -LU±+ UJ q 

n-B+np 1 + nB—riF 



Po+w 



T 



Po+uj t + OJ, 



(18) 



where ns = rifl(u g ) and n F = n F (uj±) or n F (uj T ), respectively. 

The four terms in Eq. (fT8]l correspond to the processes with the energy 
relations indicated in the denominator, i.e. the decay N — > cf>£, the production 
Ncf) — > £, the production N£ — > <fi and the production of N£(j> from the vacuum, 
as well as the four inverse reactions |17| . We are only interested in the process 
N <R- cf>£, where the decay and inverse decay are illustrated by the statistical 
factors 

1 + fiB — riF = (1 + n B ){l - n F ) + n B n F . (19) 



Our term thus reads 



1 LU. 

2u7„ ~ 



k 2 1 



n B - n F 



2m 



PO—b)±- UJ q 



(20) 
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(24) 



For carrying out the integration over the angle 77, we use 

Im = -n5(p -w±- u) q ) = -71-^5(77 - rj±), (21) 

pa - cj± - uj q kp 

where 

V± = ^ [2po^± - A/ 2 - (u£ - fc 2 ) + m 2 ] (22) 

denotes the angle for which the energy conservation po = w + holds. The 
integration over 77 then yields 

/ dr7lm(T^D*A ± ) = ^ [1 + n B (u> q ±) - n F (uj±)], (23) 

where u q ± = pa — lu± . It follows that 
r(P) = - ^- tr[(j> + M)Im E(P)] 

= 2^ Im f 2T S /(ga tr[ ( ^ + M)P^P fl ]^} 

Im I T £ / dfc dr > fc2jD * t A + (P° - + A - (Po + P'?)] } 

/" - fc 2 

-V/ dk——^—[l + n B (uj q ±)-n F (cj±)} 

x [2 Po (fc =F w±) ± M 2 ± (w 2 - fc 2 ) T mj], 

where we only integrate over regions with — 1 < 77 < 1. 

Using finite temperature cutting rules, one can also write the interaction 
rates for the two modes in a way which resembles the zero-temperature case [17) 

T±{p) = wJ a1a ^ 2 ^ 5 ^ p - k -q)\ m ^ k )\ 2 (25) 

X [1 + TIB {oj q ) - n F (ui± )] , 

where 

H 3 fc 

dfc = - d * (26) 
(27r) 3 2fc 

and dq analogously and the matrix elements are 

\M ± (P,K)\ 2 = 9 2 "\~f "± (PoTPV±) ■ (27) 

Now that we have arrived at an expression for the full HTL decay rate of a 
Yukawa fermion, we would like to compare it to the conventional approximation 
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adopted by [7J. To this end, we do the same calculation for an approximated 
fermion propagator 

1 



$ - mi 

This yields the following approximated interaction rate: 



(28) 



c(P) 



g 



32irp p 
1 

'Wo 



dk — [1 + n B (u}„) — n F (uj)][M 2 + m? — m 2 . 



dfc dg (2tt) 4 <5 4 (P - # - Q)|.M| 2 

x [1 + n B (uj q ) - njf(w)], 
where w 2 = fc 2 + ?ti 2 , uj q = pq — lo and the integration boundaries 



(29) 



fci, 



1 



p J(M 2 +mj 



2NP 

ensure — 1 < 77 < 1, where 
'/ 



) 2 - (2M TO f ) 2 T P(M^ + m| - mi) (30) 



Wp ^ 



M 



We sec that the matrix clement is 



■(M 



ml). 



This result resembles the zero temperature result 



r T = (P) 



dk - [A/ 2 



32irp p , 

with zero temperature masses me, m^. The missing factor 
1 + n B - n F = (1 + n B )(l ~ n F ) + n B n F 



(31) 



(32) 



(33) 



(34) 



accounts for the statistical distribution of the initial or final particles. As pointed 
out in more detail in [3], we have shown that the approach to treat thermal 
masses like zero temperature masses in the final state [7] is justified since it 
equals the HTL treatment with an approximate fermion propagator. However 
this approach does not equal the full HTL result. 

Concluding this calculation, a caveat has to be added: The external Ma- 
jor ana fermion will also acquire a thermal mass of order gT . Thus, if its zero 
temperature mass is smaller than that, the external fermion also needs to be de- 
scribed by leptonic quasiparticlcs to be consistent. However, in our leptogenesis 
application, the Yukawa coupling giving rise to the Majorana neutrino decay is 
much smaller than the couplings giving rise to the thermal masses of the Higgs 
boson (scalar) and the lepton (Dirac fermion) and thus the thermal mass of the 
heavy neutrino can be neglected. 

We have calculated the decay rate assuming a Majorana particle, but the re- 
sult can be very easily generalized to the case of two Dirac fermions by inserting 
the appropriate factors of two in the decay rate and the thermal masses. 
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4 Neutrino decay in Leptogenesis 



When turning to leptogenesis with 

SC = iNidrfNi - K M Ni<t>U a - ^MiNiN? + h.c, (35) 

we sum over the two components of the doublets, particles and antiparticles and 
the three lepton flavors. Thus we need to replace g 2 — s- 4(Aj,A I/ )n. Integrating 
over all neutrino momenta, the decay density in equilibrium is 

7 ^ q = / W? rNq{E)TD = i £ dEE Pf C N q r D , (36) 

where E = p n , f^(E) = {exp(E(3) — is the equilibrium distribution of the 
neutrinos and T D = [1 - f${E)} T. 

Since A„,i Q -C 1, the thermal masses are 

ml(T) = (±£ + -L<£ + \y 2 + \x) T 2 (37) 

and 

<m = [yi + L gl y. (38) 

The couplings denote the SU(2) coupling 52, the U(l) coupling gy, the top 
Yukawa coupling y t and the Higgs self coupling A, where we assume a Higgs 
mass of 115 GeV. The other Yukawa couplings can be neglected since they are 
much smaller than unity and the remaining couplings are renormalized at the 
first Matsubara mode 2ttT as explained in [?]■ 

In Fig. [5j we compare our consistent HTL calculation to the approximation 
adopted by [7] , while we add quantum statistical distribution functions to their 
calculation which then equals the approach of using an approximated lepton 
propagator l/(]j£ — me) as in Eq. (|28|) [9] . We assume the heavy neutrino masses 
to be hierarchical and evaluate the decay rate for the typical value M\ = 10 10 
GeV, which is inspired by putting M3 to the GUT scale (10 15 GeV) and assum- 
ing M1/M3 ~ 10~ 5 analogous to the quark sector. The combination of Yukawa 
couplings (Aj,Aj,)n which governs the decay rate is often parametrized by the so 
called 'effective' neutrino mass rh\ = (A„AJ,)n v 2 /Mi, where v = 174 GeV is the 
vacuum expectation value of the Higgs field. We take rh\ — 0.06 eV, inspired 
by the mass scale of the atmospheric mass splitting. However, our results can 
be generalized to all regions of parameter space. 

In the one-mode approach, the decay is forbidden when the thermal masses 
of Higgs boson and lepton become larger than the neutrino mass M < rni + m^,. 
Considering two modes, the kinematics exhibit a more interesting behavior. For 
the positive mode, the phase space is reduced due to the larger quasi-mass and at 
AI = m + (oo) + m0, the decay is only possible into leptons with small momenta, 
thus the rate drops dramatically. The decay into the negative, quasi-massless 
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Figure 5: The neutrino decay density with the one lepton mode approach 70 
and the two-mode treatment *y± for M\ = 10 10 GeV and fn\ = 0.06 cV. The 
thresholds for the two modes (+), (-) and one mode (0) are indicated. 



mode is suppressed since its residue is much smaller than the one of the positive 
mode. However, the decay is possible up to M = m^. Due to the various effects, 
the two mode rate differs from the one mode approach by more than one order 
of magnitude in the interesting temperature regime of z = T/M > 1. 

It is extremely tempting to put this result in a Boltzmann-solver and obtain 
an effect for the produced baryon asymmetry. However, in the quest for consis- 
tent treatments which capture effects of the same origin and size, other effects 
need to be included as well. At higher temperatures, when > M + rn±(k), 
the Higgs can decay into neutrino and lepton modes and this process acts as 
production mechanism for neutrinos [7]. Moreover, the CP asymmetry needs 
to be calculated taking into account the two lepton modes in order to have a 
consistent treatment. 



5 Conclusions 

As discussed in detail in [9], we have, by employing HTL resummation and finite 
temperature cutting rules, confirmed that treating thermal masses as kinematic 
masses as in [7] is a reasonable approximation. However, quantum statistical 
functions need to be included as they always appear in thermal field theory. 
Moreover, the full HTL lepton propagator shows a non-trivial two-mode behav- 
ior which is not accounted for by the conventional approach. We have calculated 
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the effect of the two modes in a general way which is applicable to any decay 
and inverse decay rates involving fermions at high temperature. Thus, this cal- 
culation is a valuable tool for other particle processes in the early universe, as 
other lcptogcncsis processes, the thermal production of gravitinos or the like. 

The behavior of the decay density of the two lepton modes can be ex- 
plained by considering the dispersion relations u>± of the modes and assigning 
momentum-dependent quasi-masses to them. The thresholds for neutrino decay 
reported in [7] are shifted and the decay density shows deviations of more than 
an order of magnitude in the interesting temperature regime T/M ~ 1. Thus 
we expect these effects to have a sizeable impact on the final baryon asymmetry. 
However, in order to arrive at a minimal consistent treatment, also the decay 
<p — > Nl at high temperatures needs to be included as well as a CP asymmetry 
that is corrected for lepton modes. In a further step, it will be interesting to 
include the effect of thermal widths in the calculations. 

As for all effects arising from thermal field theory, the effects are only im- 
portant in the weak washout regime, where lcptogcncsis takes place at high 
temperatures. We are aware of the progress that is currently being made in 
approaching the effects of quantum statistics [30H33] , quantum transport equa- 
tions |34f(5§] or other collective phenomena as e.g. the Landau-Pomeranchuk 
effect [5] ■ These efforts contribute to getting an idea of the size and impact of 
various thermal effects by approaching the extremely complex situation from 
different angles. 

Acknowledgements We would like to thank Georg Raffelt, Florian Hahn- 
Wornle, Steve Blanchet, Matthias Garny, Marco Drewes, Wilfried Buchmiillcr 
and Annika Wohner for fruitful and inspiring discussions. 

A Analytical solution for HTL lepton dispersion 
relations 

The dispersion relations of the two lepton modes are given by the poles of the 
corresponding propagator. Hence, we seek the zeros of 




( 



±fco — k ko + k 



) 



1 



±1 - 



2k ko — k 



(39) 



The equations D± = can be transformed by the substitutions 



fco + k 

ko — k 
ko — k 



1 



(40) 



(41) 



k + k 
k 2 



■>■+ 



(42) 



c := — . 
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This yields 

D± = ±- — - — (-2c-l + x±-\nx±). (43) 
c x± — 1 

Further introducing 

s := - exp(-2c - 1) (44) 

leads to 

D± = 1 [x ± + ln(- S ) - l na;± ] . (45) 

1 + m(— s) x± — 1 

Since the prefactor does not have poles for the values of K we are looking at, 
solving D± = amounts to solving 

x± + ln(-s) - In x± = 0, (46) 

which in turn means 

s = -x±c-' x± . (47) 

This is the defining equation of the Lambert W function [40|41j . thus the 
solution reads 

x± = -W(s). (48) 



According to the definition in Eq. (|44[) 

- 1/e < s < 0, (49) 

thus the two real branches of the Lambert function, Wo and W—i, correspond 
to the two solutions we seek. In the range given by Eq. (|49| Wo > — 1 and 
W-i < — 1. For fco > k we have x+ > 1 and x_ < 1. Hence, the physical 
solutions for x± read 

x + = -W- 1 (s) and x-=-W Q (s). (50) 

The corresponding results for u± are then given by 

u+ = k w-M + i (51) 



W (s) - 1 



(52) 



Wo(s) + 1 
Making use of the relations [42] 

W ^ 1 (z)+hi(Wo,-i(z))=lnz, (53) 
one can directly prove the result by plugging Eqs. (jBTj) and ([52^1 into Eq. (|3"9"|) . 
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